Many body effects on c- axis properties: out of plane coherence and bilayer splitting. 
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The out-of-plane hopping in layered metals with strong electronic correlation is analyzed theo- 
retically. By studying the effects of in-plane interactions on the interlayer tunnelling we investigate 
one of the oldest unresolved problems of high-Tc cuprates, the so-called bilayer splitting. 

The strong electron-electron interactions within each layer reduce the quasiparticle weight, and 
modify the hopping between layers. We analyze the effect of the in-plane correlations on the inter- 
layer hopping using a Renormalization Group scheme already applied to the problem of interlayer 
coherence . The bilayer band splitting acquires a significant temperature and doping dependence, 
and can be completely suppressed when the Fermi energy coincides with a Van Hove singularity. 



I. INTRODUCTION 

The strong anisotropy and anomalous c- axis prop- 
erties of the high-temperature cuprate superconductors 
are among the most relevant features of these layered 
materials and have greatly conditioned the theoretical 
scenario^. Although intensively investigated, the out-of- 
plane behaviour of the cuprates is not yet well under- 
stood. The role of the multiple layers and its relevance 
in the critical temperature Tc and optimal doping values 
of different compounds, is not settled yet^. 

One of the most interesting issues has been the in- 
teraction between Cu02 planes in materials with two 
planes per unit cell. As predicted by early band struc- 
ture calculations, the interaction between planes would 
split their electronic structure in bonding and antibond- 
ing bands, the so-called bilayer band splitting. This split- 
ting was expected to be maximum at (tt, 0) and van- 
ish along (0,0) — (7r,7r), by symmetry effects^. The ab- 
sence of the BS on early angle-resolved photoemission ex- 
periments on Bi2Sr2CaCu20s+s (Bi22f 2) aroused great 
expectationi, in conection with the interlayer pair tunnel- 
ing modeli. Recently, angle-resolved photoemission spec- 
troscopy (ARPES) data on bilayer cuprate Bi2212 have 
reported a bilayer band splitting (BS) between bond- 
ing and antibonding bands in the Brillouin zone (BZ) 
antinodal region of overdoped samples, above and below 
Tjj^. The bilayer energy splitting, due to the intrabi- 
layer coupling, is of about f 00 meV in the (tt, 0) region 
in the normal state and of about 20 meV in the super- 
conducting state. The momentum and energy behavior 
of the BS agree qualitatively with the bilayer Hubbard 
model calculations^. The report of BS as well in under- 
doped and optimally doped Bi22f 2 samples, not resolved 
before, by ARPES measurements with different photon 
energies and high momentum and energy resolution*, 
should be seriously considered and put constrains in ex- 
isting theories and data analysis2ii2iii. The existence of 
BS in all doping ranges seems to confront with the co- 
herence issue. In this context, a crossover in the phase 
diagram between the low temperature, overdoped side 



with coherent electronic excitations, and the high tem- 
perature, underdoped side where this coherence is lost, 
has been reported from ARPES data and resistivity mea- 
surements of Bi22I2-'^^, in good agreement with theoreti- 
cal predictions^'^, and with optical conductivity resulta^^ 
from which the coherent-to-incohcrcnt crossover as a 
function of doping has been reported. In the overdoped 
region, conventional Fermi liquid behavior is expected 
while in the underdoped region correlation effects would 
deviate from the conventional physics yielding exotic be- 
havior, in agreement with the most common doping- 
temperature phase diagram for the cuprates. A crossover 
from insulating-like, at high temperatures, to metallic- 
like character at low temperatures in the direction per- 
pendicular to the layers, while being metallic within the 
layers, has been recently reported by Valla et al from 
ARPES measurements in layered, non-superconducting 
materialsi^, the change of effective dimensionality being 
correlated with the existence or non-existence of coher- 
ence within the layers. The splitting of the CM02-plane 
states in the whole doping range would imply the coher- 
ence of electronic states in the c axis at least to the intra- 
cell distance between the two planes. At present, consen- 
sus on the existence of BS has been reached only in the 
overdoped regime, while in optimally doped and under- 
doped regime there is not a clear agreement. Recentlyi^ 
high-resolution ARPES experiments on Bi22I2 covering 
the entire doping range, and with fotons of 22 eV and 
47 eV (ARPES matrix element effects for fotons of 47eV 
strongly enhance the antibonding component of the bi- 
layer in the (tt, 0) region of the BZ) suggest that the 
size of the intracell coupling is only weakly dependent on 
the doping level, therefore some coherent coupling exists 
between the Cu02 planes of the unit cell even in un- 
derdoped system. To explain the c-axis transport mea- 
surements where an insulating pc behavior is reported in 
underdoped samples, ref^^ argue that intracell perpen- 
dicular hopping will dominate the BS while the intcrcell 
perpendicular hopping will control the c-axis resistivity 

Pc 

It is the purpose of this paper to analyze the role played 
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by many-body effects on the interlayer coupling of high- 
Tc superconductors and its effects on the bilayer splitting. 



II. MODEL AND METHOD 

In spite of the great theoretical effort in high- Tc su- 
perconductivity, there is not yet a theory of consen- 
sus. The Hubbard hamiltonian is among the most stud- 
ied models in the hight-Tc superconductors field, since 
it explains the main low-energy physics results of the 
Cu02 planes, which are the common blocks to all fam- 
ilies and where it is assumed to lie the main physics of 
the cuprates. The insulating parent compound is ob- 
tained at half-filling with antiferromagnetic correlations, 
while upon doping the known phase diagram is obtained. 
The stripe phase as well as d-wave superconductivity are 
well described within the Hubbard modeUi. Further- 
more, in the overdoped regime, where Fermi-liquid-like 
behavior seems to hold, ARPES data of Bi2212 qual- 
itatively agree with bilayer Hubbard model calculation 
results for the energy splitting-», a experimental intrabi- 



fluctuations of the system: 



layer hopping t^,. 



44 meV was obtained. 



However, the anisotropy of the normal state transport 
properties: electron motion in the c-direction is inco- 
herent in contrast with the metallic behavior of the in- 
plane electrons as probed by the different and pab 
resistivitiesiSiiS, cannot be understood in the current the- 
oretical framework. 

The anomalous out of plane behavior of the cuprates 
has led to the suggestion that conventional Fermi liquid 
theory fails in these compounds^. Both theoreticali^fli^ 
and experimental^^ work has remarked the relevance of 
the properties of the direction perpendicular to the Cu02 
planes. 

An alternative explanation of the emergence of inco- 
herent behavior in the out of plane direction has been 
proposed in terms of the coupling of the interlayer elec- 
tronic motion to charge excitations of the systemJ^. This 
approach implicitly assumes that electron electron inter- 
actions modify the in-plane electron propagators in a non 
trivial way. 

We will show that, even in the clean limit, many body 
effects can suppress the coherent contribution to the out 
of plane electron hopping. The clean limit is defined as 
that in which the length scale, L, over which electrons 
remain coherent within the layers diverges. 

The simplest formulation of the method replaces the 
excitations of the system (such as electron-hole pairs) 
by a bath of harmonic oscillators with the same excita- 
tion spectrum. This approach can be justified rigorously 
in one dimension, and is always an accurate description 
of the response of the system when the coupling of the 
quasiparticles to each individual excitation is weakS^. 

In the following, we will assume a local interaction be- 
tween electrons close to the Fermi level, and the charge 



H-int — C - 



k 



(1) 



where Ci creates an electron at site «, and describes the 
charge fluctuations of the environment, which are to be 
described as a set of harmonic modes. The Hamiltonian 
of the system is approximated as: 



TLe-b — 'Helec + Ti-env + 'Hint (2) 



where Heiec describes the individual quasiparticles, 'Henv 
stands for the set of harmonic oscillators which describe 
the environment, and Hint defines the (linear) coupling 
between the two. The b\ are boson creation operators, 
the tij describe the electronic hopping processes, and the 
information about the interaction between the electron in 
state i and the system is defined by the the function— 

Using second order perturbation theory and eq.([3lj we 
can writeSiSi: 



J.(u;) = ^y/(k)Imx(k,o.) 



(3) 



where x(k, lu) is the Fourier transform of the density- 
density response of the system, (pg(i)/5_jj(0)). The inter- 
action in eq.(Q is spin independent. Other, more compli- 
cated, couplings can also be taken into account, provided 
that the appropriate response function is used. 

The influence of the electron-boson coupling on the 
electron propagators can be calculated to all orders if 
the state i is localized, that is, neglecting the hopping 
terms in eq.©. We find: 



exp < — did 1 — e 



Ji{u)) 



(4) 



where {c\{t)ci{t'))Q ~ g^ei{t-t ) -g ^-^c Green's function in 
the absence of the interaction. The method that we use 
assumes that eq. also holds in a system with extended 
states. For a standard metallic system, we must insert 
{cl{t)ci{t'))o ~ l/(t - t') in eq.©. It can be shown that 
this approximation is exact at short times, W <ti {t — 
t')~^ <C A, where W is an energy scale related to the 
dynamics of the electrons, and A is the upper cutoff in 
the spectrum of the environment. 

The time dependence in eq.(@} is determined by 
lini(^^o x(k, w). In a gapless, metallic system, we have 
x(k, a;) ^ a{]<i)\Lj\. This behavior, when inserted in 
eq.Q, leads to: 



{c\{t)c,{t')) 



(t 



(5) 
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where: 



dkV^Ck)a(k) 



(6) 



k|<L- 



where L is the scale of the region where the tunneling 
process takes place. The value of L is hmited by the 
length over which the phase of the electronic wavefunc- 
tions within the layers is well defined. We assume that, in 
a translationally invariant system, there is no dependence 
on the position of the local orbital, i. This result implies 
that the frequency dependence of the Green's function, 
in a continuum description, can be written as: 



lim G(f — r', cx 



(7) 



We can now use eq.© to analyze the interlayer tunnel- 
ing by applying Renormalization Group methods. The 
simplest case where this procedure has been used is 
for the problem of an electron tunneling between two 
states, i and j, which has been intensively studiedSS'2&. 
We integrate out the high energy bosons, with energies 
A — dA < cjfc < A and rescaled hopping terms are de- 
fined. As mentioned earlier, eq.© is valid for this range 
of energies. The renormalization of the hoppings is such 
that the properties of the effective Hamiltonian at en- 
ergies Lo A remain invariant. If the hoppings tij are 
small, any physical quantity which depends on them can 
be expanded, using time dependent perturbation theory, 
in powers of: 

tl{cl{t)c,it)clit'Mt')} « t|.(4(t)c,(t'))(c,(t)4(t')) 

(8) 

The integration of the high energy modes implies that the 
terms in eq.® are restricted tot < A~^, or, alternatively, 
the time unit have to be rescaled^l, r' — re'^^^^, where 
T ~ A~^. Using eq.jSJ, the condition of keeping the 
perturbation expansion in powers of the terms in eq.® 
invariant implies that: 



dA 

t e~ 



2 ^i^(2+2Q) 



(9) 



which can also be used to define the scaling dimension of 
the hopping terms. Finally, 



9(t,,/A) 



(10) 



dl A 

where I = log(Ao/A), and Aq is the initial value of the 
cutoff. 

This approach has been successfully used to 
describe inelastic tunneling in different situations 

The analysis which leads to eq. Hl()|l can be generalized 
to study hopping between extended states, provided that 
we can estimate the long time behavior of the Green's 
function, G'(k,< - t') = (ct(i)c£ (*'))■ This function is 
related to the local Green's function, which is written in 
eq.©, by: 



^ hm G(f - F, w) = / d'^kG(k, u) 

|r— r' |-^0 



(11) 



where D is the spatial dimension. In the cases discussed 
below, the interaction is instantaneous in time, and the 
non interacting Green's function can be written as: 



Go(k,c.)oc-^^ 



(12) 



where 2 = 1,2. In the following, we assume that the 
interacting Green's function has the same scaling prop- 
erties, with the factor uj~^ replaced by u}~^ in ea. p2|l . 
where 5 depends on the interactions. This can be shown 
to be correct in perturbation theory to all orders, in the 
models studied below, because the corrections depend 
logarithmically on lo (it is a well known fact for the Lut- 
tinger liquid). Then, using eqs. 0, ((TT|l and ifT^ . we 
obtain: 



G(k,cj) oc ^ 



(13) 



and J-{u) is finite. Thus, from the knowledge of the real 
space Green's function, using eq.Q, we obtain a, which, 
in turn, determines the exponent a-\- D / z which charac- 
terizes G(k,Ci;). Generically, we can write: 

Gi^Lo) ^ \uj\'^-^ (14) 

where the subindices I, e stand for localized and extended 
wavefunctions. In terms of these exponents, we can gen- 
eralize cg. dlOII to tunneling between general states to: 



9(t:;/A) 

dl 



' A 



(15) 



Before proceeding to calculations of &i and (5e for various 
models, it is interesting to note that, in general, the re- 
sponse function of an electron gas in dimension _D > 1 
behaves as lim^^g x(k, a;) ~ |a;|/|k|, so that, from 

eq.®, limL—oo « ^ L'-^^^\ Thus, for D > I, the contri- 
bution of the inelastic processes to the renormalization 
of the tunneling vanishes for delocalized states, L — » oo. 



A. Van Hove singularities in the density of states. 

The dispersion relation of two-dimensional electronic 
systems with a square lattice, i.e. the Cu02 planes of 
the high- Tc superconductors, usually present Van Hove 
singularities. The Fermi surface of most hole-doped 
cuprates is close to a Van Hove singularity. The possible 
relevance of this fact to the superconducting transition 
as well as to the anomalous behavior of the normal state 
was put forward in the early times of the cuprates and 
gave rise to the so called Van Hove scenario^. We will 
assume that the metallic layers are well described by elec- 
trons in a square lattice, and that the Fermi level is close 
to the (tt, 0){A) and (0, 7r)(i?) points of the Brillouin Zone 
(BZ). Close to these points, the dispersion relation can 
be parametrized as: 



eA,B(k) 



2m 



2m. 



(16) 
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where rrix and my are parameters which can be estimated 
from the band structure of the model. In the foUowing, 
we win consider the renormahzation of the interlayer tun- 
nchng associated to these regions in the BZ. 

The response function at low energies and small 
wavevectors has been computed in^: 



Imx(k, uj) 



E 

i=A,B 



1 



47re,(k) 



(17) 

where ei(k) is the dispersion relation 116(1 . 

The long time dependence of the Green's func- 
tion is determined by the low energy behavior of x- 

limoj^olm x(k,a;) ^ Hi=A,B\^\/^ii^)- We assume 
that the interaction between the electrons and the den- 
sity fluctuations is short ranged as before. The diver- 
gence of Imx when eA,B = implies that the inte- 
gral in eq.© diverges logarithmically as L — > oo, as 
in the two previous cases, irrespective of the details 
of the interaction, V^(k). Because of this divergence, 
it is convenient to shift slightly the chemical potential 
£f away from the saddle poinlj^, evh- A finite value 
of |eF — Svh\ implies the existence of a length scale, 



£vh\ 



Lq ^ [Max(rn3;, mj,)|ei? 

the k integrals in eq.||HJl. 

Using a local potential, we find 



1/2 



which regularizes 



a^iU/W)^\og\L/Lo) 



(18) 



where W is an energy scale of the order of the width of 
the conduction band. The dependence of a on L goes as 
log^(L), as in other physical quantities in this model24. 

In the two dimensional model, z = 2, so that di = a, as 
estimated above, and — a — 1. The divergence of a im- 
plies that tunneling between localized and also between 
extended states is suppressed at low temperatures. In ad- 
dition, the effective electron-electron coupling, U, grows 
at low energies or temperatures, until a scale at which the 
system is unstable and a phase transition takes place—. 
This effect enhances the suppression of interlayer hop- 
ping. 



III. BILAYER BAND SPLITTING 

To study de bilayer band splitting (BS) we apply the 
method described above and consider the interlayer tun- 
nelling between two Cu02 planes within the Van Hove 
scenario. We will consider the electronic states near the 
two inequivalent saddle points of the square lattice. The 
lowest order corrections to the bare couplings included in 
the t-t'- Hubbard Hamiltonian, will give the anisotropic 
screening of the interactions^, by the Kohn-Luttinger 
mechanism, in contrast to the isotropic three-dimensional 
metals. The parameters of the model will be the inter- 
layer hopping, the distance of the Fermi level £j7 to the 
Van Hove singularity evh given by the band structure, 
the band width W and the strength of the interaction a. 



A{o}) 



w 



CO 



FIG. 1: Schematic representation of spectral function A(lj) 
versus frequency for the strong coupling case. The band width 
W indicates the high-energy cutoff. 



The effective Hamiltonian will be defined at each scale of 
the couplings. The \£vh — £f\ defines the low energy cut- 
off while the bandwidth W defines the high-energy cutoff 
of the Renormahzation Group approach. The hopping 

renormahzation gives the te// « t ^ ^^^^^^^ ^ 

We analyze the phase diagram relating the strength 
of the couplings to the doping level. The underdoped 
regime would correspond to the strongly correlated sys- 
tem. The interacition a > 1, the Van Hove singular- 
ity lies at the Fermi level so that the low-energy cutoff 
\evH ^ £f\ = 0; no BS occurs at this regime, a broad 
structureless background will extend in the scale of W as 
it is shown in Fig 1. 

The system is instable at low energies, and low tem- 
perature phase transition to the superconducting phase 
are highly probable. The underdoped-optimally doped 
regime would correspond as well to strongly correlated 
system, a > 1, but in the frequency region below the low- 
energy cutoff a renormalized bilayer splitting will appear, 
since quasiparticles would be renormalized. The renor- 
mahzation factor depends on the cutoffs as 



Z = 



\ w 



Between the low and high energy cutoffs, a broad 
background without any structure would be obtained, as 
schematically represented in Figure 2. 

The overdoped regime, as stated above, would corre- 
spond to the weakly correlated case, a < 1. Electron cor- 
relations are screened and, as schematically represented 
in Figure 3, the renormahzation of the bilayer splitting 
would be weaker than in the preceding case. An incoher- 
ent tail of the spectral weight is shown with a power law 
decay at high energies, ^(i^) ~ sin{Tra/2) always 
below the high-energy cutoff. 
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IV. CONCLUSIONS. 




FIG. 2: Same as Fig.l with renormalized bilayer splitting 
at frequencies below the low-energy cutoff, \evH — £f\- The 
band width W indicates the high-energy cutoff. 
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FIG. 3: Same as Fig.l with bilayer splitting at frequencies 
below the low-energy cutoff, \evH ~ £f\- The spectral weight 
decays at high energies. The band width W indicates the 
high- energy cutoff. 



The effects of electron-electron correlations on the out- 
of-plane properties of layered systems have been stud- 
ied. We have discussed the suppression of interlayer tun- 
neling by inelastic processes in two dimensional systems 
in the clean limit. Our results suggest that, when per- 
turbation theory for the in-plane interactions leads to 
logarithmic divergences, the out of plane tunneling ac- 
quires a non trivial energy dependence. The conductance 
goes to zero as T ^ if the Fermi level of the interact- 
ing electrons lies at a Van Hove singularity. Coherence 
would be suppressed and the bilayer splitting would van- 
ish as a consequence of the strong correlations. When 
\tVH — 7^ 0, a strongly renormalized bilayer split- 
ting occurs below \evH — ^f\ which marks the low-energy 
cutoff of the system. These two situations correspond 
to the underdoped and optimally doped regions in the 
cuprate scenario. In the overdoped case, quasiparticles 
would be weakly renormalized and the bilayer splitting 
will be present. 

Thus, we have shown that insulating behavior in the 
out of plane direction is not incompatible with gapless or 
even superconducting in-plane properties, although the 
in plane properties are also markedly different from those 
of an ordinary Fermi liquid. 
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